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1 Introduction 



In previous work pQ, two of us determined non- perturbative D°, D + , D*+, D+, and 
A+ fragmentation functions (FFs), both at leading order (LO) and next-to- leading or- 
der (NLO) in the modified minimal-subtraction (MS) factorization scheme, by fitting the 
fractional-energy spectra of these hadrons measured by the OPAL Collaboration [215] in 
e + e _ annihilation on the Z-boson resonance at the CERN Large Electron-Positron Col- 
lider (LEP1). Apart from untagged cross sections, they also measured the contributions 
arising from Z — > bb decays. This enabled the authors of Ref. [1] to obtain specific FFs 
for the transitions c,b — > D°, D + , D* + , D+, A+. The strategy adopted in Ref. p] was very 
similar to the one underlying Ref. [I], in which also ALEPH data [S] were fitted, and 
Ref. [B] . The FFs obtained in Ref. pQ were used as input for a NLO study |7J of charmed- 
meson hadropro duct ion in pp collisions, which yielded reasonable agreement with data 
collected by the CDF Collaboration in run II at the Tevatron [8]. 

Recently, new data on charmed-meson production with much higher accuracy have been 
presented by the Belle Collaboration [9] at the KEK Asymmetric Electron-Positron Col- 
lider for B Physics (KEKB) and the CLEO Collaboration [TO] at the Cornell Electron- 
Positron Storage Ring (CESR). These data offer us the possibility to determine the non- 
perturbative initial conditions of the FFs much more accurately. Furthermore, the large 
span in center-of-mass (cm.) energy (\fs) ranging from 10.5 GeV [9pTJ] way up to 
91.2 GeV [2,3,5] provides us with a powerful lever arm to test the Dokshitzer-Gribov- 
Lipatov-Altarelli-Parisi (DGLAP) [11] evolution of the FFs. These new FFs will enable 
us to improve our theoretical predictions for the charmed-meson hadroproduction cross 
sections [71TT2] to be compared with the CDF data [8]. 

The data from Belle and CLEO are located much closer to the thresholds \/s = 2m c and 
y/s = 2m& of the transitions c — > H c and b —> Hb, where H c and Hf, stand for generic c or b 
hadrons, respectively, than those from ALEPH and OPAL. It might thus be a questionable 
approximation to treat the partonic cross sections for e + e~ — > c + X and e + e~ —> b + X 
in the massless approximation, with m c = rrib = 0, as was done in Refs. [HUE], where 
LEP1 data were fitted. Therefore, we take into account the finite quark mass corrections 
of the form m 2 /s (m = m c ,mb) in the partonic cross sections to test their significance. 
Similar studies based on perturbative FFs [13] may be found in Refs. [T^ffTS] . 

The outline of this paper is as follows. In Section [21 we describe the theoretical formalism 
of single-hadron inclusive production in e + e~ annihilation. After reviewing the massless 
case in Section 12.11 we explain how to include the full mass corrections in Section I2.2L 
Lengthy expressions are relegated to Appendices [A] and [Bl Specifically, we list the elec- 
troweak quark charges in Appendix [A] and the mass-dependent coefficient functions in 
Appendix [B] In Section 12.31 we consider the electromagnetic initial-state radiation (ISR) 
that is inherent to the Belle and CLEO data and explain how to efficiently accommodate 
the ISR corrections in our fits. In Section [HI we present several alternative FF sets. They 
are obtained from global fits to Belle % CLEO [E], ALEPH [5], and OPAL (2EJ data, 
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and from separate fits to the B-factory (Belle plus CLEO) and Z-factory (ALEPH plus 
OPAL) data, both for m^O and for m = 0. Finally, in Section HJ we present a summary 
and our conclusions. 



2 Formalism 



We study the inclusive production of a single charmed hadron H c , with mass m#, in e + e 
annihilation via a virtual photon (7) or Z boson, 

e+ + e- ^( 1 ,Z)^H C + X, (1) 

where X stands for the residual final state, which goes unobserved. Specifically, we 
concentrate on the cases H c = D°, D + , D* + . In the following, we explain how to calculate 
the cross section of process ([1]) at NLO in the parton model of QCD, both in the zero- 
mass (ZM) approach, where all quark masses are neglected, and in the general-mass (GM) 
approach, where the c and b quarks are taken to be massive. We denote the four-momenta 
of the virtual gauge boson and the H c hadron by q and pn, respectively, so that s = q 2 
and m 2 H = p%, and introduce the scaling variable x = 2(p# • q)/q 2 - We call the energy of 
H c and the angle of its three-momentum w.r.t. the beam axis in the cm. frame E and 9, 
respectively. Then, x = 2E/y/s measures the energy of H c in units of the beam energy. 

For unpolarized beams and observed hadrons, the cross section of process ([I]) at a given 
value of ^/s can only depend on E and 9. Since the virtual boson has spin one, the most 
general form of the differential cross section then reads 

d 2 a Z M 2m da T 3 . 2n da L 3 n da A , n , 
-(l + cos 2 6»)- h-sm 2 6>- h-cos#-— . (2) 



dxdcos^ 8 dx 4 dx 4 dx 

The three terms on the right-hand side are the transverse, longitudinal, and asymmetric 
contributions, respectively. The first two are associated with the corresponding polar- 
ization states of the virtual boson with respect to the direction of the observed hadron. 
The asymmetric contribution is due to the parity-violating interference terms and is not 
present in QED. The transverse and longitudinal parts are normalized so that 

da 7 d 2 a da T da L 

dx J dxd cos 9 dx dx 

-1 



In the parton model, each component dcr p /dx (P = T,L) on the right-hand side of 
Eq. ([3]) can be written, up to power corrections, as a sum of convolutions of partonic cross 
sections derf (y, fi, [if)/ dy, where a = g, u, u, . . . , b, b is the fragmenting parton with four- 
momentum p a , y = 2(p a ■ q)/q 2 , and \i and [if are the renormalization and factorization 
scales, respectively, with FFs D a (z,fif), where z = x/y is the fraction of energy passed 
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on from parton a to hadron H c in the cm. frame, as 



dc , \ f dy dcrf . . „ / x 



dx J y dy \y 



£ -ZsttoM t )DA-,,i t , (4) 



where the values of y m i n and y max are subject to mass effects to be discussed below. At 
NLO, /!/ defines the scale, where the divergence associated with collinear gluon radiation 
off a massless primary quark or antiquark is to be subtracted. 



2.1 ZM approach 



At NLO in the MS scheme, the cross sections of the relevant partonic subprocesses are 
given by 



( Nl/ )AW ) = iV> iyl + Al) \ 5(1 -y) + nA,l) 



dy 

dOg 

dy 



2tt 



"f 



a s (n) 



2tt 



Here, iV c = 3 is the number of quark colors; 



(5) 



00 



47TQ; 5 

~37" 



(6) 



with a being Sommerfeld's fine-structure constant, is the total cross section of e + e~ — > 
fi + n~ for massless muons; V qi and A q . are the effective vector and axial-vector couplings 
of quark q { to the photon and the Z boson including propagator adjustments, which are 
listed in Appendix [A] Pf^ 1 are the LO timelike splitting functions [TT] . 



Pt T Xv) = c F 



3 An n , 1 +y 2 

J 

y 



and the coefficient functions read 

9 2 



C,(y) = Ct 



+ 4 



In y 



1-3/ 



2 3 y v y; 2 Vi-?/ 

(l + y)[21ny + ln(l-y)] 



Mi - y) 

1-3/ 



1 J_ fl _y)2 

C,(y) = ^[2 lny + ln(l - y)], 



(7) 



where Cp = (N% — 1)/(2N C ) = 4/3 and the plus distributions are defined as usual. We 
evaluate the strong-coupling constant a s using the two-loop formula with rif = 5 quark 
flavors. We identify fj, — [if — y/s , so that in Eq. (jSJ) the terms proportional to ln(s//ij) 
vanish. We observe that C g (y) < for any value of y, so that the gluon contributes 
destructively to Eq. (jHJ). The bounds of integration in Eq. (TJJ are y m i n = x, y mSuX = 1, and 
we have yfpu < x < 1, where p# = Am 2 H / s. 

For later use, we also list the total hadronic cross section at NLO, 

nj 
8=1 

The x dependences of the FFs are not yet calculable from first principles. However, once 
they are given at some initial fragmentation scale /io, their fif evolution is determined by 
the DGLAP evolution equations [TT] . 

i 

^D a (x, n) = ?M'£j ^ (y , a s ({i))D b (i, „/) . (10) 

Specifically, we use /io = rn c for a — g,u, u, d, d, s, s, c, c and po = m b for a — b,b. Our 
task is thus to construct a model for the z dependences of D a (z, /i ), which upon evolution 
to fif = y/s fit the data at that cm. energy. 

The above formalism is identical to the one that is routinely used in the literature for 
the inclusive production of single light hadrons [T7]. The non-zero values of the c- and 
6-quark masses only enter through the initial conditions of the FFs, and the mass of the 
heavy hadron sets the lower bound on the scaling variable x. 

2.2 GM approach 

We derived the partonic cross sections at NLO for non-zero quark masses adopting the 
on-shell definition of the latter and found agreement with Ref. [18J. We take the pole 
masses of the c and b quarks to be m c = 1.5 GeV and m& = 5.0 GeV, respectively. The 
finite-mass corrections are generally of order m 2 / s (m = m c , nit,). They can be sizeable for 
sufficiently small values of y/s. In our study, where the smallest value of yfs is 10.52 GeV, 
they reach at least 2% for the c quark, but might be much larger depending on the 
coefficient of m 2 /s in the LO cross section. On the other hand, they are quite substantial 
for the b quark at values of y/s just above the bb production threshold. 

In the following, we consider finite-mass effects of order m 2 /s only in the production dy- 
namics, but not in the decays of the produced H c mesons, which are affected by kinematic 
power corrections. This would be beyond the scope of our analysis and is not required 
because the experimental analyses actually provide production cross sections. 



i + ?Mc P * 

2tt 2 



(9) 
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In order to expose the connection with the ZM limit and to establish the finite subtraction 
terms, which are needed for the evaluation of the GM cross section in the MS factorization 
scheme, we report the relevant cross section formulae for the GM case. We largely adopt 
the notation of Ref . [18] . 



At NLO, there is a contribution from real gluon radiation. The gluon can either be 
included in the hadronic system X or it can act as the leading parton initiating a hadron 
jet that contains the H c meson. Therefore, we also need to consider at the parton level 
single-gluon inclusive production with X including a cc or bb pair. 

The partonic cross sections of single-heavy-quark inclusive production read 



dy 



N c a 



{v \y,p) + A 2 q F^\y,p) 



where p = Am 2 / s with m = m c , mj. Notice that the vector and axial-vector contributions 



differ by finite-mass terms. At NLO, the coefficient functions F P (y,p), with P 
and u = v, a, in Eq. ( TTTT) may be decomposed as [IS] : 



F { p\y, P ) 



6(l-y)B$>(p) + 



2vr 



6(l-y)s£\p) 



i-y 



T,L 



(12) 



The LO terms read 



B?>( P )=f3, B^\ P ) = P 4, 



4 a) (p) = /? 3 



4 a) (p) = o, 



(13) 



where (3 = y/1 — p. In the massless limit, we have B^\p) = B^\p) = 1 and B^\p) = 
B£\p) = 0. The NLO terms, S { p\p) and R { p\y,p), may be found in Appendix [B] The 
partonic cross section of single-gluon inclusive production reads 



?(<*)/ 



N c a 



2vr 



V q ^\y,p)+AlG^(y,p) 



(14) 



with the coefficient functions G <y p\y, p) given in Appendix! 



Deviating from the notation 



of Ref. 



we included the factor Cp in the functions Sp (p), Rp (y,p), and Gp(y,p) 



to simplify the comparison with the ZM approximation. Through this comparison, we 
recover the so-called perturbative FFs [13], from which the subtraction terms for the 
conversion from the NLO calculation with finite quark masses in the on-shell scheme [IS] 
to the GM variable-flavor-number scheme [71fT2] are constructed. The result thus obtained 
contains all the finite-mass terms and, at the same time, smoothly approaches the ZM 
result in the limit m — > 0. 

Thanks to our specific choice of hadronic and partonic scaling variables, x and y, respec- 
tively, Eq. (jlj) applies to the GM case as it stands. The bounds of integration in y and the 
allowed x range now depend on the partonic subprocess and on the fragmenting parton 
a. In the case of heavy-quark fragmentation, we have y m i n = max(x, yfp), 2/max = 1, and 
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\fp~H < x < 1. If a gluon fragments and there is also a heavy-quark pair in the final state, 
then we have y min = x, y max = /3 2 , and \fp~H < x < (3 2 . If there are only massless partons 
in the final state, then we have y m i n = x, y max = 1, and \Jp~n < s < 1 as in Section [2TTT . 

The experimental data at y/s = m z , collected by OPAL [2J3] and ALEPH [5J at LEP1, 
come in the form da/dx as a function of the scaling variable x introduced above. As 
explained above, the maximum x range is y/pH < x < 1. On the other hand, Belle 
[9] and CLEO pH] present their data as distributions da/dx p in the scaled momentum 
x p = p/Pmax = \/ (x 2 — Ph) /(I — Ph) , with allowed values < x p < 1. The conversion 
formula reads 

= (15) 



with x = J {I - p H )x 2 p + Ph- 

We conclude this section with a remark concerning parton model kinematics of frag- 
mentation in the presence of finite quark and hadron masses. In the picture where the 
fragmenting parton a creates a jet that includes the observed hadron H, the parton vir- 
tually p 2 must exceed the hadron mass square m 2 H . In fact, if the scaling variable z is 
defined in terms of light-cone momenta, as z = p H /Pa = (p°h +Pn)/(Pa + Pa) [19J , we have 
Pa > vr? H j z [20]. Since the outgoing partons are taken to be on-shell in the parton model, 
this inequality becomes m 2 > m 2 H /z. In our case, this is only satisfied for the transition 
b —>■ H c in the constrained z range m 2 H / m 2 < z < 1 . 

At this point, we find find it instructive to insert a digression on the massive kinematics 
of fragmentation, a topic which has received very little attention in the literature. Let us 
view the fragmentation process b — > H c as a decay b — > H c + X , where the invariant mass 
nix of the hadronic system X, which we may treat as one effective particle, can be tuned 
in the range < nix < m — ran- Starting from the rest frame of the decay, we perform 
a Lorentz boost along the three-momentum of H c , such that the energy of H c becomes 
P°h = x y/s/2j where \/s is the e + e~ cm. energy. This corresponds to the situation where 
H c is emitted collinearly from b, carrying energy p® = y\/s/2, in the laboratory frame. 
For given values of x and mj, we have 



2m 2 H 



(m 2 + m 2 H — m 2 x) x — \J A (m 2 , m H , m 2 x ) (x 2 — pu) 



(16) 



where A (a, b, c) = a 2 + b 2 + c 2 — 2(ab + be + ca). Owing to the constraint y/p < y < 1 on 
y and the one on mj specified above, x must lie in the range J~pu < x < x max , with 



■'■max 2 



2m 2 



[m 2 + m 2 H + (3 (m 2 - m 2 H )] . (17) 
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For a given value of x, the allowed y range is thus y m i n (x) < y < y max (x), where 

[jp if J~pH < x < 

ymm{X) — <; m 2 + m 2 , 



where 



Up 



2m 2 H 



min(yo, 1) H y/pH < x < (m + m H )^/^ 
min(x^,lj if (m + m H )^/^ < x < x max 



(m 2 + m 2 ^) x — (m 2 — m 2 H ) ^/x 2 — pn ■ (19) 



For ^Js > (m 2 + m 2 H )/mH, the scaling variable z = x/y appearing in the b — > H c FF lies 
in the range 2m 2 H /(m 2 + m 2 ^) < z < x max . 

If m# <C m is a good approximation, we may simplify these expressions by putting 
m H = 0. Then, we have < x < x max , with x max = (1 + (3)/2, and 

^/p if < x < & 



!hmn[ ' r> ^x + f x iff <x<x 



max 



ymax(x) = 1. (20) 

Consequently, we have < z < x max . Finally, if also m = 0, we recover the ZM situation, 
< x < 1 and x < y < 1, so that < 2; < 1. 

In Fig. (U the kinematically accessible region in the (x, y) plain is shown for our frag- 
mentation model of b — > assuming LEP1 experimental conditions, with -^i = m^, 
m = 5.0 GeV, and m# = 2.01 GeV (solid line). It does not exhaust the region en- 
compassed by yfpu < x < 1 and max(x, ^J~p) < y < 1 that we adopt for the eval- 
uation of Eq. (Tj0) in the case of b — > H c (dashed line). There is a wedge missing at 
small values of x and large values of y, i.e. small values of z. The minimum value of z, 
2m 2 H /(m 2 + m 2 H ) Pd 0.28, is reached at the minimum value of x. However, in the bulk 
of the missing wedge, we have z < mnjm « 0.40. As we shall see in Sectional this 
excluded z range actually accommodates the peak of the b — > H c FF (see also Table [HD- 
Moreover, our simple fragmentation picture cannot describe the transitions a —>■ H c for 
a = u,d,s,c,g. For these reasons, we abandon it at this point. A crucial conceptual 
drawback of the parton model applied to the inclusive production of heavy hadrons is 
that fragmenting partons are taken to be on-shell. We escape this problem by imagining 
that the hadronic system initiated by the fragmenting parton a not only receives color, 
but also energy and momentum from the rest of the event. In this way, the accessible 
region in the (x, y) plane is expanded to become the one underlying our analysis. By the 
same token, the masses of the fragmenting quark and the observed hadron can be treated 
independently of each other, the former appearing in the partonic cross sections and the 
latter in the hadronic phase space factor. As explained above, we take m = < m# in 
the ZM approach and m, ran > in the GM approach. In Section El we also study the 
case m = ran = for comparison. 
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2.3 Electromagnetic initial-state radiation 



The cross sections of inclusive single-hadron production in e + e annihilation measured 
by Belle 0, CLEO [10J, ALEPH 0, and OPAL [2 : 3J naturally include electroweak cor- 
rections, which were not subtracted in the data analyses. The bulk of these corrections 
is due to the effect of electromagnetic radiation emitted from the incoming electrons and 
positrons. This ISR is suppressed by a factor of a, but enhanced by the large logarithm 
log(s/rn.g), where m e is the electron mass. At Belle and CLEO energies, the hadronic 
cross section decreases with increasing invariant mass of the hadronic system. Since ISR 
reduces the hadronic mass, it leads to an increase in cross section. The shape of the FF is 
also changed, since a fraction of the events takes place at lower hadronic invariant mass. 
The impact of ISR on the determination of FFs has already been considered in Ref. [Hj 
and has been found to be non-negligible for the analysis of the Belle and CLEO data. 

The most straightforward way to correct for ISR would be to incorporate into the general 
expression for the e + e~ — > H c + X cross section without photon radiation, presented in 
Section 12.11 the corrections due to photon radiation off the initial-state leptons and use 
the resulting expression to fit the FFs to the Belle and CLEO data. This procedure would 
involve several additional numerical integrations in each iteration of the fitting procedure 
and thus dramatically slow down the latter. In the following, we explicitly derive an 
approximation formula, which reduces the number of integrations to a manageable level 
and is still rather precise. Our procedure differs from the one used in Ref. [T3j, where the 
ISR corrections were subtracted in an iterative way from the experimental data before 
the actual fit. Our procedure is numerically more involved, but offers the advantage that, 
at at the end, the ISR corrections precisely refer to the final x distribution resulting from 
the fit. 

The dominant ISR corrections are conveniently incorporated using the structure-function 
approach, in which the photon emission is taken to be collinear to the incoming e ± leptons 
[2"Tf2"2"] . In analogy to the factorization formula of the collinear parton model, the ISR- 
corrected differential cross section (1(Jisr(p+, P-), where p± denote the four-momenta of 
the incoming e leptons, is obtained by convoluting the uncorrected differential cross 
section da(p + ,p-) with radiator functions D e ±{x±, s), one for each incoming leg, which 
measure the probabilities for the e leptons to retain the fractions x± of their energies 
after the emission of ISR, as 



dcrisR(jo+,P-) = / dx+ / dx- D e +(x + , s)D 



[x. 



s)da(x + p + , X-p. 



(21) 



Using the method by Gribov and Lipatov |23j, the leading logarithms can be resummed 
to all orders, leading to the expression for D e ±(x, s) in Eq. (7) of Ref. [22J. The structure- 
function approach was mostly applied to total cross sections in the literature. In this case, 
one integration can be carried out independently of the considered process and leads to 
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the luminosity function 

1 1 

H e + e -(r,s) = J dx + J dx_<5(r — x + x_)/} e +(x + , s)D e -(x_, s) (22) 
o o 

in t = x + x_. The Gribov-Lipatov-resummed expression for it may be found in Eqs. (8) 
and (9) of Ref. [22], in terms of the variable x = 1 — r - However, if the cross section is 
differential w.r.t. variables, whose differentials are not invariant under boosts along the 
beam axis, as in our case, the situation is more involved [21]. In this case, we obtain from 
Eq. (121 p the following master formula: 



dcriSR 
dx 



[x, s) = J dx + dx_ dx' dcos^' S(x — x(x + , x_, x', cos9'))D e +(x + , s)D e - (x_, s) 



dV 



x 



dx' d cos 9 



-(x , cos9', x + X-s), 



(23) 



where the primed variables refer to the hadronic cm. frame, which is reached from the 
e + e~ cm. frame through a Lorentz boost with velocity (3 = (x + — x_)/(x + + x_), and 



nn I nn ff <~Y> 
•JU I JU _|_ j JU ^ JU 



', cos 9') = 7 Lfttf + fi^rx' 2 - p H cos & 



(24) 



with j — 1/yl — /3 2 E] Integrating Eq. fTSBl over cos#', we have 



daisR , 
dx 



[x,s 



dx' 



a+ / (a;'— a— /x+) 



a+/(x'—a—/x—) \ 



/ 1 

J dx + J dx_ + J dx_ / dx 

\x/x' a— /(x'— a+/x+) x/x' a—/(x'—a-\-/x—) J 



(25) 



where 



2D e +(x+, s)D e - (x-, s) d 2 a ( , 2x — (x< + x)x' 

x v + ' — — x, v + . ' =, r.s 

| x+ - x_ | a/x /2 - pn/r dx'd cos 9' y (x + - X-)^x 12 - p H /r 

x ± a/x 2 - p H 



a± 



(26) 



Notice that the second pair of integrations enclosed within the parentheses in Eq. ( 1251) 
merely duplicates the first one if D e +(x, s) = D e -(x, s). 

Detailed numerical inspection of our specific application reveals that, to very good ap- 
proximation, we may substitute in Eq. (|2"3"|) 

D e +(x+,s)D e -(x-,s) « 5(1 - x + )H e+e -(x_,s), (27) 



1 We denote the relativistic boost velocity and its 7 factor by /? and 7, respectively, because (3 and 7 
are reserved for other quantities in this paper. 
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which trivially satisfies Eq. (I2"2"j) . so as to eliminate the second term within the parentheses 
in Eq. (I25p and to save one integration in the first one. This implies that the whole ISR 
is emitted by the electron alone, while the positron stays idle, or vice versa. A further 
simplification may be obtained by casting Eq. (1251) into a form that contains on the r.h.s. 
(da/dx')(x' , s), i.e. the cross section at the e + e~ cm. energy integrated over the polar 
angle. To this end, we use the approximations 

ix' C0S9', TS) « -(1 + COS 2 9')^-(x' TS) 



dx' d cos 9' 8 da;' 



3 /-, 2 &(ts) da 



The approximation in the first line of Eq. ( 1281) may be justified by observing that, in 
Eq. (|2J), dcr L /dx is suppressed for m 2 <C s, being zero at LO in the ZM approach, and that 
the contribution from da A /dx vanishes upon integration over cos 9'. The approximation 
in the second line of Eq. (|28p faithfully describes the leading power-like dependence on the 
cm. energy, but disregards the scaling violations of the FFs, which are just logarithmic. 
For simplicity, we evaluate the factor <j(ts)/(j(s) in Eq. ( 1281) at LO in the GM approach 
neglecting the contribution from Z-boson exchange, as 



a(s 




a{rs) 1 l + p/(2r 



(29) 



On the one hand, the NLO corrections largely cancel out in this cross section ratio; on 
the other hand, the Z-boson contribution is suppressed for a/s <C mz, as explained in 
Appendix [A] Substituting Eqs. (1271) and (1281) into Eq. (1251) . we obtain our working formula 

*?™(xs)k [dx' — (x's) dr 2H ^-^ S ^ h 1 + co * ^ 

dx M-Jte^x^) J dr (i _^^____ 8 (l + cos ^ )(t(s) , (30) 

x a-/(x'— a + ) 

where 

a, 2x - (1 + t)x' f . 

cos 9 = = 31 

{l-r)^x' 2 -p H /r 

and a(rs)/a(s) is given in Eq. ( |29l) . In Eq. ( |30|) . we insert the NLO expression for 
(da/dx')(x',s). Strictly speaking, we should then also evaluate Eq. ( |29i) at NLO. How- 
ever, the omitted correction is insignificant, as explained below Eq. ( 1291) . We employ 
this formalism also in the ZM approach, except that we then set p = in Eq. (I29p for 
consistency. 

The effect of ISR on the x p distribution of e + e~ — > D + + X measured by Belle [9] and 
CLEO [IU] is studied for the GM approach in Fig. [21 where the result of the ISR-corrected 
fit to these data is compared with the corresponding result where the ISR corrections are 
subtracted. As already observed in Ref. [Hj, we find that the spectrum is shifted to larger 
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values of x p , i.e. it becomes harder and lower at the peak when the ISR corrections are 
subtracted. This is expected, since ISR reduces the available hadronic cm. energy, which 
softens the spectrum and increases the cross section. 

At LEP1 energy, the shift of the spectrum due to ISR is negligible. This may be under- 
stood by observing that ISR shifts the cm. energy available for the hard scattering to 
values below the Z-boson resonance, where the cross section is greatly reduced. In the 
analysis described below, we thus only include ISR corrections in our theoretical descrip- 
tion of the Belle [H] and CLEO [TU] data, while we neglect them in connection with the 
ALEPH and OPAL [2B] data. 



3 Results 

As experimental input for our fits, we use the x p distributions of D°, D + , and D* + 
production in the continuum at a/s = 10.52 GeV for 0.08 < x p < 0.94 from Belle [9] 
and for 0.20 < x p < 0.95 from CLEO |10j . which we correct for ISR as described in 
Section [231 and the x distributions of D°, D + [2], and D* + [3"]l5] production on the Z- 
boson resonance at \fs = 91.2 GeV from ALEPH and OPAL. We received the Belle data 
in numerical form via private communication [25]. Belle [9j also provide data from the 
T(5S) resonance outside the 5-meson decay region, which we leave aside. 

In Refs. [2f3f5] . the absolute cross section distributions in x are normalized to the total 
hadronic cross section and include the branching fractions of the decays used to identify 
the H c mesons, namely D° — > K~tt + , D + — > K~tt + tt + , and D* + — > D°tt + followed by 
D° — > K~7i + , respectively. Therefore, we multiply our predictions by l/a tot , where, for 
simplicity, we evaluate a tot from Eq. © exploiting the insignificance of quark mass effects 
at LEP1 energy, and divide the experimental data by the respective decay branching frac- 
tions. For consistency, we adopt the very values of the latter that are used in Refs. [2|3ll5] . 
These read B(D° -»• K~x + ) = (3.84±0.13)% and B(D + -> K-it + it + ) = (9.1 ±0.6)% [26] 
for Ref. 0, B(D*+ -» D°n+) = (68.3 ±1.4)% and B(D° -> R-tt + ) = (3.83 ±0.12)% [27] 
for Ref. |3], and B(D*+ -> D°n + ) = (68.3 ±1.4)% and B(D° -> K~ti + ) = (3.85 ±0.09)% 
[28] for Ref. [5], respectively. The 6-tagged samples of Refs. [2|3] are treated in the 
same way as the full cross sections, while the one in Ref. [5] needs to be multiplied by 
Rbfip — > D^), where Rb = T(Z — > bb)/T(Z — > hadrons) is the fraction of 6-tagged events 
in the full hadronic sample and f(b — > D*^) is the probability of a b quark to hadronise 
into a D*^ 1 meson. We adopt the value Rbf{b — > D* + ) = (4.66 ± 0.51)% determined in 



We take a in Eq. © to be the running fine-structure constant, which is particularly 
important because it appears there in squared form. At y/s = 10.52 GeV, we have 1/a ~ 



132 [29]. Of course, this effect cancels out in the normalized cross sections considered in 




As already mentioned in Section [1} we perform for each hadron species a combined fit to 



Ref. [5]. 




Refs. [2113115] . 
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Table 1: Values of fit parameters for the D° meson resulting from the Belle/CLEO, OPAL, 
and global fits in the GM approach together with the values of x 2 achieved. 





Belle/CLEO-GM 


OPAL-GM 


global- GM 


N c 


1.51 • 10 7 


4.42 • 10 4 


8.80 • 10 6 


a c 


1.56 


1.52 


1.54 


7c 


3.64 


2.83 


3.58 


N b 


13.5 


13.5 


78.5 


a b 


3.98 


3.98 


5.76 


lb 


0.921 


0.921 


1.14 


X 2 


3.15 


0.794 


4.03 



the Belle, CLEO, ALEPH, and OPAL data (global fit), a separate fit to the 5-factory data 
(Belle/CLEO), and a separate fit to the Z-factory data (ALEPH/OPAL). We perform 
fits both for m ^ (GM) and for m = (ZM). 

We parameterize the z distributions of the c and b quark FFs at their starting scales /zo 
as suggested by Bowler [30], as 

D^(z,n ) = Nz~^ 2 \l - z)^l\ (32) 

with three free parameters, N, a, and 7. This parameterization yielded the best fit 
to the Belle data [9] in a comparative analysis using the Monte-Carlo event generator 
JETSET/PYTHIA. 

Specifically, our fitting procedure is as follows. At the scale /i/ = m c = 1.5 GeV, the 
c-quark FF is taken to be of the form specified in Eq. ( |32i) . while the FFs of the light 
quarks q (q = u, d, s) and the gluon are set to zero. Then these FFs are evolved to 
higher scales using the DGLAP equations in Eq. (1101) at NLO with n/ = 4 active quark 

flavors and an appropriate value A^jL of the asymptotic scale parameter. When the scale 

JVl o 

reaches the threshold value fif = nib = 5.0 GeV, the bottom flavor is activated and its 
FF is introduced in the Bowler form of Eq. (I3"2l . The evolution to higher scales is then 
performed with n f = 5 and the value A^ is properly matched to A^L Including A^ 
among the fit parameters, it turns out to be feebly constrained by the fit. Therefore, we 
adopt the value A^ = 221 MeV from Ref. [31] and adjust the value of A^ accordingly, 

to be A^L = 321 MeV. 

We first describe the Belle/CLEO fits in the GM approach. In the Belle and CLEO data, 
all charmed hadrons coming from 5-meson decays are excluded, so that there is no need 
to include b — > H c fragmentation. On the other hand, the ALEPH and OPAL data each 
come as two sets: the sample of H c hadrons produced by the decays of b hadrons from 
Z — > bb (b tagged) and the total sample of H c hadrons, also including those from direct 
production in Z — > cc and from light-quark and gluon fragmentation (total). Since we 
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Table 2: Values of fit parameters for the D + meson resulting from the Belle/CLEO, 
OPAL, and global fits in the GM approach together with the values of x 2 achieved. 





Belle/CLEO-GM 


OPAL-GM 


global-GM 


N c 

a c 

7c 


5.66 • 10 5 
1.15 
3.39 


2.82 • 10 4 
1.49 
2.92 


5.67- 10 5 
1.16 
3.39 


N b 

a b 

lb 


18.8 
4.71 
1.17 


18.8 
4.71 
1.17 


185 
7.08 
1.42 


X 2 


1.30 


0.509 


1.99 



Table 3: Values of fit parameters for the D* + meson resulting from the Belle/CLEO, 
ALEPH/OPAL, and global fits in the GM approach together with the values of x 2 
achieved. 





Belle/CLEO-GM 


ALEPH/OPAL-GM 


global-GM 




1.33 ■ 10 7 


4.58 • 10 4 


1.10 • 10 7 


a c 


0.992 


1.38 


1.07 


1c 


3.84 


3.00 


3.81 


N b 


6.67 


6.67 


14.0 


a b 


3.28 


3.28 


3.85 


lb 


1.04 


1.04 


1.14 


X 2 


3.74 


2.06 


6.90 
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wish to test the Belle/CLEO FFs fits through comparisons to LEP1 data, we must include 
the b — > H c transitions in an appropriate way. For this purpose, we first fit the D° and 
D + data from OPAL and the D* + data from ALEPH and OPAL. The resulting values of 
the fit parameters and of x 2 P er degree of freedom, x 2 , for the D°, D + , and D* + mesons 
are given in Tables HHSl respectively. The goodness of these fits may also be judged from 
Figs. E](a)-(c), respectively. We observe from Fig. [3^c) that the D* + data from ALEPH 
and OPAL are only moderately compatible, as was already noticed in Ref. [I]. This 
explains why the corresponding fit has a larger value of x 2 than those for the D° and D + 
mesons. 

In a second step, we use the values of N b , a b , and % thus obtained as rigid input for 
the fits to the Belle and CLEO data. While the b — *> H c transitions are excluded from 
the final states as explained above, their FFs still influence the fits through the DGLAP 
evolution from [if = m h to /i/ = *J~s. These fits yield new values for N c , a c , and j c , which 
are also included in Tables HH3] together with the values of x 2 achieved. The x 2 values 
of the fits discussed so far are all acceptable, except perhaps for the fits to the D° and 
D* + data from Belle and CLEO, which yield x 2 values in excess of 3. In Figs. IH(a)-(c), 
the D°, D + , and D* + data from Belle and CLEO are compared with the respective fit 
results in order to assess the goodness of the latter. We observe from Figs. Eta)-(c) that 
the agreement is best in the D + case, which is also reflected in Tables HH3J 

We now turn to our global fit, which uses all available D°, D + , and D* + data, from Belle, 
CLEO, ALEPH, and OPAL. The resulting values of the fit parameters and of x 2 are a l so 
included in Tables HH31 The D°, D + , and D* + data are compared with the respective 
theoretical results based on the global fit in Figs. [5H3, respectively. Comparing Figs. [5H7] 
with Figs. [3] and HI we find the most striking difference to be that the global fit describes 
the total H c samples from ALEPH and OPAL, especially the D° and D + samples from 
OPAL, less well in the large-x range, for x^>0.6, than the ALEPH/OPAL fits. This may 
be understood by observing that, in the global fit, the c — > H c FFs, which dominantly 
contribute in the large- x range, are mostly constrained by the more precise Belle and 
CLEO data and are appreciably increased in that x range compared to their counterparts 
from the ALEPH/OPAL fits. On the other hand, the modifications of the b — > H c FFs are 
less significant and do not worsen the agreement with the ALEPH and OPAL data in any 
visible way. The observation that the c — > D* + FF fitted to Belle and CLEO data does 
not yield a satisfactory description of the ALEPH data was also made in Ref. [H] , where 
it was speculatively linked to the presence of large non-perturbative power corrections. 
New experimental data at intermediate energies, between those of the B and Z factories, 
which are not expected to become available in the foreseeable future, would certainly shed 
more light on this potential anomaly. 

In order to study the impact of finite c and b quark masses on our fits, we repeat them in 
the ZM approach, where m c = m b = 0, except in the definition of the starting scale /i . We 
still have m# ^ 0, so that the x distributions have finite lower endpoints and differ from 
the corresponding x p distributions in shape. The resulting values of the fit parameters and 
of x 2 are listed for the D°, D + , and D* + mesons in Tables HHEJ respectively. Comparing 
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Table 4: Values of fit parameters for the D° meson resulting from the Belle/CLEO, OPAL, 
and global fits in the ZM approach together with the values of x 2 achieved. 





Belle/CLEO-ZM 


OPAL-ZM 


global-ZM 




1.03- 10 7 


3.43 • 10 4 


1.04- 10 7 


a c 


1.48 


1.48 


1.50 


1c 


3.60 


2.80 


3.60 


N b 


13.4 


13.4 


80.8 


a b 


3.96 


3.96 


5.77 


lb 


0.923 


0.923 


1.15 


X 2 


3.25 


0.789 


4.66 



Table 5: Values of fit parameters for the D + meson resulting from the Belle/CLEO, 
OPAL, and global fits in the ZM approach together with the values of x 2 achieved. 





Belle/CLEO-ZM 


OPAL-ZM 


global-ZM 


a c 

1c 


7.30 • 10 5 
1.12 
3.43 


2.62 • 10 4 
1.48 
2.91 


7.31 • 10 5 
1.13 
3.43 


N b 
a b 

lb 


19.0 
4.71 
1.17 


19.0 
4.71 
1.17 


163 
6.93 
1.40 


X 2 


1.37 


0.507 


2.21 



Table 6: Values of fit parameters for the D* + meson resulting from the Belle/CLEO, 
ALEPH/OPAL, and global fits in the ZM approach together with the values of x 2 
achieved. 





Belle/CLEO-ZM 


ALEPH/OPAL-ZM 


global-ZM 


N c 


1.05 ■ 10 7 


2.80 • 10 4 


1.14 • 10 7 


a c 


0.929 


1.33 


1.03 


1c 


3.82 


2.93 


3.82 


N b 


6.52 


6.52 


14.9 


a b 


3.25 


3.25 


3.87 


lb 


1.04 


1.04 


1.16 


X 2 


3.69 


2.04 


7.64 



16 



Tables HHH1 with Tables HH3J we observe that the inclusion of finite quark masses reduces 
the x 2 values of the global fits by 11-16% and also tends to reduce the x 2 values of the 
Belle/CLEO fits, except for the case of D* + mesons, where the difference is insignificant. 
As expected, the quality of the ALEPH/OPAL fits are practically unaffected by finite- 
quark- mass effects, which provides a retrospective justification for the use of the ZM 
approach in Refs. [HUE]- 

From the comparison of the fit parameters in the GM and ZM approaches presented in 
Tables HH3] and 0HH1 respectively, it is hard to judge by how much the FFs of the two 
approaches actually differ as functions of z at a given value of fif. As an example, we 
thus display in Figs. [Sfa) and (b) the z dependences at [if — 10.52 GeV of the c — ► D + 
and g — > D + FFs, respectively, of the Belle/CLEO fits in the GM and ZM approaches. 
We notice that there is only little difference between the GM and ZM results. This 
observation may be understood from Fig. [9j where the ZM result for the x p distribution 
of e + e~ — > D + + X at \fs = 10.52 GeV evaluated with the GM FFs is compared with the 
proper GM result. In fact, the finite-m c correction to the hard-scattering cross section only 
amounts to a few percent. A similar observation was made in Ref. [32] using perturbative 
FFs P3]. 

In the above implementation of the ZM approach, mu is identified with its physical values. 
It is interesting to study the impact of the finite-m# correction. To this end, we repeat 
the Belle/CLEO-ZM fit for the D + meson putting also m# = 0, which implies that x p = x 
and (da/dx p )(x p ) = (da/dx)(x), as may be gleaned from Eq. ffT5"]) . In order to obtain 
an acceptable value of x 2 i we exclude the six data points with x p < 0.2 from the fit. 
Furthermore, we require that r > p D + in Eq. ( l23l) . in which pn = is put otherwise, 
to ensure that the hadronic energy after ISR is above the production threshold. The 
resulting c — > D + and g — > D + FFs are also shown in Figs. E](a) and (b), respectively. We 
notice that the former significantly differs from its counterpart in the proper ZM approach, 
its peak being reduced in size and shifted to a lower value of z, while the latter is only 
moderately affected. This modification of the FFs is compensated by a reciprocal change 
in the line shape of the x p distribution, as may be seen from Fig. [91 which also contains 
the result of the ZM approach with mj? = evaluated with the Belle/CLEO-GM FFs. In 
fact, the peak position and height are substantially increased relative to the evaluations 
with rnn 7^ 0. The situation is similar for the D° and D* + mesons. 

Besides the c —>■ H c and b —>■ H c FFs themselves, also their first two moments are of 
phenomenological interest. They correspond to the branching fractions, 



B Q (f Jlf )= J dzD Q (z,fif), (33) 



max( v /piF,,z C ut) 



where Q = c,b, and the average fraction of energy that the H c meson receives from the 
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Table 7: Values of c — > if c and 6 — > if c branching fractions at /i/ = 2m fe , 10.52 GeV, and 



FF set 


H c 


5 C (10.52 GeV) 


B c (m z ) 


B b (2m b ) 


B b (m z ) 


Belle/CLEO-GM 


D° 


0.525 


0.611 


0.146 


0.492 




D+ 


0.232 


0.269 


0.0590 


0.168 




D*+ 


0.211 


0.249 


0.0696 


0.206 


ALEPH/OPAL-GM 


D° 


0.493 


0.591 


0.146 


0.491 




D+ 


0.185 


0.220 


0.0590 


0.167 




D*+ 


0.200 


0.247 


0.0695 


0.206 


global-GM 


D° 


0.522 


0.608 


0.140 


0.490 




D+ 


0.230 


0.268 


0.0512 


0.157 




D*+ 


0.206 


0.245 


0.0716 


0.212 


Belle/CLEO-ZM 


D° 


0.534 


0.622 


0.146 


0.490 




D+ 


0.235 


0.273 


0.0592 


0.167 




D*+ 


0.215 


0.254 


0.0695 


0.205 


ALEPH/OPAL-ZM 


D° 


0.489 


0.587 


0.146 


0.489 




D+ 


0.185 


0.221 


0.0591 


0.166 




D*+ 


0.201 


0.248 


0.0694 


0.204 


global-ZM 


D° 


0.527 


0.614 


0.141 


0.488 




D+ 


0.234 


0.272 


0.0517 


0.157 




D*+ 


0.209 


0.248 


0.0718 


0.210 
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Q quark, 

1 

(z)cj{Vf) = Bgifif) j zDq(z, [if), (34) 

max( v ^7^",2 cut ) 

where the cut z cnt = 0.1 excludes the problematic z range where our formalism is not valid. 
As may be seen from Figs. [3] and HI there are no experimental data at z < z cut either. 
Tables [7J and M contain the values of Bq([if) and (z)q([if), respectively, for Q = c,b and 
H c = D°, D+, D*+ at fi f = 2m b , 10.52 GeV, and m z for the Belle/CLEO, ALEPH/OPAL, 
and global fits in the GM and ZM approaches. We observe from Table [7] that, within the 
GM and ZM approaches, the values of B c (fif) from the ALEPH/OPAL fits are somewhat 
smaller than those from the Belle/CLEO fits, by less than 10% for the D° and D* + 
mesons and by approximately 20% for the D + meson. The corresponding results from 
the global fits tend to lie between those from the Belle/CLEO and ALEPH/OPAL fits, 
but closer to the former. The GM and ZM approaches yield very similar results. As 
for Bb(fif), the differences between three fits and two approaches are minor. We note 
that the values of B^m^) have to taken with a grain of salt because, in conrast to the 
c — > H c FFs, the b — ► H c FFs are not directly constrained by low-energy data. Looking 
at Table El we see that the values of {z)q(fif) are shifted towards smaller values through 
the DGLAP evolution in [if, as expected, and that quark-mass effects are insignificant 
here. As for (z) c ([if), the values from the ALEPH/OPAL fits fall 5-9% below those from 
the Belle/CLEO fits, which are only slightly larger than those from the global fits. As 
for the ALEPH/OPAL and global fits, this trend may already be noticed by comparing 
Figs. [3(a), (b), and (c) with Figs. 0(b), EJb), and[7J(b), respectively. On the other hand, 
the differences between the various fits are marginal for (z)b([if)- 

Of course, the Q — > H c FFs and their moments depend on scale, scheme, order, and 
implementation issues such as the functional form of the ansatz at the starting scale [to 
and the value of [Iq itself, and thus do not represent physical observables by themselves. 
Nevertheless, comparisons of the quantities Bq([if) and (z)q([if) defined in Eqs. ([33]) and 
f[3"4"|) . respectively, with their experimental counterparts determined from the measured x 
distributions are instructive, as they reveal in how far the Q —> H c transitions actually 
dominate the cross section distributions. Moreover, these quantities stringently charac- 
teristise the lineshape in x of the Q — ► H c FFs at a given value of [if and simplify the 
comparisons with our previous FF sets [41(6] and those to be introduced by other authors. 

The values of the branching fractions and average energy fractions of the Q — > H c tran- 
sitions measured by Belle and CLEO at ^Js = 10.52 GeV and by ALEPH and OPAL at 
y/s = m z are collected in Tables M and [TUl respectively. For comparison, we present in 
Table [TT1 the counterparts of B c {m z ), B h {m z ), (.z) c (10.52 GeV), and (z) c (m z ) extracted 
from the cross section distributions based on the global fits in the GM approach, which 
are shown in Figs. EHTJ 

In the remainder of this section, we compare our favourable FFs, from the global fit in the 
GM approach, with those from Refs. [Tf4~!6"], which were determined through fits to ALEPH 
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Table 8: Values of average energy fractions for c — > H c and b — > H c transitions at 
fif = 2m b , 10.52 GeV, and m z . 



FF set 


H c 


(z) c (10.52 GeV) 


(z) c (m z ) 


(z) b (2m b ) 


(z) b {m z ) 


Belle/CLEO-GM 


D° 


0.623 


0.479 


0.470 


0.273 




D + 


0.629 


0.484 


0.470 


0.293 




D* + 


0.659 


0.503 


0.508 


0.305 


ALEPH/OPAL-GM 


D° 


0.591 


0.450 


0.470 


0.273 




D + 


0.596 


0.455 


0.470 


0.293 




D* + 


0.614 


0.462 


0.508 


0.305 


global-GM 


D° 


0.621 


0.477 


0.453 


0.274 




D+ 


0.629 


0.484 


0.451 


0.288 




D* + 


0.655 


0.499 


0.501 


0.306 


Belle/CLEO-ZM 


D° 


0.624 


0.480 


0.471 


0.274 




D+ 


0.632 


0.486 


0.470 


0.293 




D*+ 


0.661 


0.504 


0.509 


0.306 


ALEPH/OPAL-ZM 


D° 


0.591 


0.450 


0.471 


0.274 




D+ 


0.596 


0.455 


0.470 


0.294 




D*+ 


0.613 


0.461 


0.509 


0.306 


global-ZM 


D° 


0.623 


0.479 


0.454 


0.275 




D+ 


0.631 


0.486 


0.452 


0.289 




D*+ 


0.657 


0.500 


0.501 


0.308 



Table 9: Values of c — > H c and b —> H c branching fractions extracted by ALEPH [5] and 
OPAL [2|3] at y/s = m z from their measured cross section distributions. 



Q - 




ALEPH 


OPAL 


c — 


> D° 


0.559 ±0.022 


0.605 ± 0.040 


c — > 


D+ 


0.238 ±0.024 


0.235 ±0.032 


c — > 


D* + 


0.233 ±0.015 


0.222 ± 0.020 


b^ 


D *+ 




0.173 ±0.020 
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Table 10: Values of average energy fractions for c — > H c transitions extracted by Belle 
[9] and CLEO [ID] at ^/s = 10.52 GeV and by ALEPH [5] and OPAL [2][3] at ^/s = 
m z from their measured cross section distributions. The values for Belle and CLEO 
are obtained by converting the corresponding average momentum fractions quoted in 
Refs. pnU] . respectively. 



H c 


Belle 


CLEO 


ALEPH 


OPAL 


D° 


0.640 ±0.002 


0.640 ± 0.005 




0.487 ±0.014 


D+ 


0.647 ±0.001 


0.650 ±0.007 




0.483 ±0.019 


D*+ 


0.682 ±0.001 


0.682 ±0.006 


0.488 ± 0.008 


0.515 ±0.009 



Table 11: Counterparts of B c {m z ), B b (m z ), (;z) c (10.52 GeV), and (z) c (m z ) extracted 
from the cross section distributions based on the global fits in the GM approach. 





ll B c (m z y 


ll B b (m z y 


"(z) c (10.52 GeV)" 


u {z)c(m z )" 


D° 


0.628 


0.515 


0.632 


0.509 


D+ 


0.276 


0.164 


0.640 


0.516 


D*+ 


0.252 


0.221 


0.666 


0.532 



[S] and OPAL [215] data in the ZM approach parameterizing the c — » H c and b — > H c 
FFs using the Peterson [33J and Kartvelishvili-Likhoded [31] ansaetze, respectively. In 
Ref. p], the initial scale for the DGLAP evolution was taken to be fio = m c ,mb as in 
the present paper, while it was chosen as fiQ = 2m c ,2mb in Refs. [4|6] . As in Figs. E^a) 
and (b), we consider the c — > D + and g — > D + FFs at fif = 10.52 GeV. The comparison 
is presented in Figs. fTOT a) and (b). From Fig. [TUT a). we observe that our global-GM 
c — > D + FF significantly differs from those of Refs. [TfB"] in lineshape, which essentially 
reflects the strong pull of the Belle and CLEO data and the difference between the Bowler 
and Peterson parameterisations. From Fig. [TOT b). we see that our global-GM g — > D + FF 
is similar to the one of Ref. [1] , while the one of Ref. [6] is strongly suppressed, especially 
in the lower z range. As explained in Refs. [TUT] , this may be attributed to the reduced 
length of the evolution path. The situation is similar for the D° and D* + mesons. 



4 Conclusions 

Previous determinations of non-perturbative charmed-hadron FFs in the parton model 
of QCD [TPfB] were based on data from the Z-boson resonance, so that the effects of 
finite quark and hadron masses were greatly suppressed and could safely be neglected. 
The advent of precise data from the B factories offers us the opportunity to further 
constrain the charmed-hadron FFs and to test their scaling violations. However, this 
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motivates the incorporation of quark and hadron mass effects, which are then likely to be 
no longer negligible, into the formalism. The GM variable-flavor-number scheme, which 
we previously advocated [71TT2] . provides a rigorous theoretical framework for this and is 
employed here for the first time to determine FFs of heavy hadrons. 

Specifically, we determined here new FFs for D°, D + , and D* + mesons through global fits 
to all available e + e~ annihilation data, from Belle [9], CLEO [ID], ALEPH [5J, and OPAL 
[2|3] . In contrast to the situation at the Z-boson resonance, the x distribution of the cross 
section for continuum production is appreciably distorted by the effects of electromagnetic 
ISR, which we, therefore, took into account. For comparison, we also performed fits to 
individual data sets. We found that the global fits somewhat suffer from the fact that the 
Belle and CLEO data tend to drive the average x value of the c — > H c FFs to larger values, 
which leads to a worse description of the ALEPH and OPAL data. Since the b — > H c 
FFs are only indirectly constrained by the Belle and CLEO data, their form is only feebly 
affected by the inclusion of these data in the fits. In order to assess the significance of 
finite-mass effects, we repeated the fits in the ZM variable-flavor-number scheme. As 
expected, the inclusion of finite-mass effects tends to improve the overall description of 
the data, by reducing the x 2 values achieved. Specifically, hadron mass effects turned out 
to be more important than quark mass effects. In fact, they are indispensable to usefully 
describe the low-x p tails of the measured cross sections. 

A FORTRAN subroutine that evaluates the FFs presented here for given values of z and \if 
may be obtained from the authors upon request via electronic mail. 



Acknowledgment 

The work of T.K., B.A.K., and G.K. was supported in part by the German Federal 
Ministry for Education and Research BMBF through Grant No. 05 HT6GUA and by the 
German Research Foundation DFG through Grant No. KN 365/7-1. The work of T.K. 
was also supported in part by the DFG through Graduate School No. GRK 602 Future 
Developments in Particle Physics. 

A Effective electroweak charges 

The effective electroweak charges appearing in Eqs. ([5]), ([9]), ( |TT1) . and f fl^l) are given by 

K = e l e \ + 2e e t> e e^pi(» + (v 2 e + a 2 e ) u 2 p 2 (s), 

Al = (v 2 e +a 2 e )alp 2 (s), (35) 

where Vf = (T 3 / — 2e/ sin 2 9 w )/(2 sin 9 W cos 9 W ) and aj = T 3 f /(2sin9 w cos9 w ) are the 
vector and axial- vector couplings of fermion /, with fractional electric charge e/ and third 
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component T 3 f of weak isospin, to the Z boson, and 



s(s - m|J 



pl[S) ( s -m%) 2 + m%r 



2 ' 

Z) T" ""Z 1 z 
2 



P2(s) = 5 - -I) 2 + m|r| (36) 

are propagator functions. Here, 9 W is the weak mixing angle and Y z is the total decay 
width of the Z boson. For small energies, \fs <C mz, the propagator functions pi(s) and 
p 2 (s) are negligible. 



B Single heavy-quark inclusive cross sections at 0(cx s ) 

In this appendix, we list the NLO coefficient functions appearing in Eqs. ( [121) and ( fT4l) . 
We cast our results in a form similar to Ref. [18] , except that our formulas include the 
overall factor Cf- We start by introducing the short-hand notation: 



T x = l- X, 



7%. 



£(z,p) = In 



p — 2x — 2^/ x 2 — p 
p-2x + 2^x 2 - p 



(37) 



The coefficient functions due to virtual-soft corrections to inclusive single heavy-quark 
production read: 



SP (p) = p 2 S$ } (p) + 2C F p0 z In 



S<£\p) = 0. 



c f 



(2-P) 



4 In - In — ^4 - 4 Li 2 
p 1-/3 



+ -vr 2 + In' 



1+/3 
1-/3 



+ Li 2 



4/3 



(l-/*) s 



+ (10 - 8p) In - 4/3 - 8/3 In - 

x p p 

2^ ln r^' 



?2 c (f ) 



1 + /3 
1-/3' 



1-/3 
2/3 

- Li 2 



2 In" 



4/3 



2/3 



/3 



(38) 
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The coefficient functions due to real corrections to inclusive single heavy-quark production 
read: 



R^\x,p)=C F 



aA 2 - P 



P{ J+ {Ar x + pf + T *V + T *>{ 1 ATx + p 



+ 



Ri\x,p) = C F 



+ 



p 2 (2 - r 2 ) + p(2x 3 - 7x 2 - 1) + 2x 2 {l + x 2 ) 
2(x 2 - p) 

-p(l-p)-r,(r,-2p) + 2T f (1 + T - 
p 3 + p 2 (4r ;c -3)+p(3a; 2 -l) 



2 (a; 2 - p) 



<(x,p) > , 



B$\x,p)=C F 



4 



3 2 rt 2r 2 (l + t x Y 



+ 



V 7 ^ 2 - P 
r 2 (l + r :c )(5r x + 4) 
4r x + p 

-2p 3 + 8p 2 x + px 2 (2r x - 9) - p + 2x 2 (l + x 2 ) 



Rt\x,p)=C F 



2t 2 



2(x 2 - p) 

8r x (l + r x f 2(1 + 7", 



+ 



\/x 2 - p 
r x p(p + r 2 - 2) 



p + x-5- 



(At x + p) 2 At x + p 



£(x,p)j , 

t 2 - 8r x - 2) 



(39) 



2 (:r 2 - p) 

The coefficient functions due to real corrections to inclusive single gluon production read: 
G ( p(x,p)=C F {2 



l + (l-x) 2 l-x p 2 
'- + P- 



X 



x 2x 



1-/9, 



- 4 



l-x 



x 



P% 

X 



Gi\z,p)=C F 
G ( *\x,p)=C F {2 



-^lni±^ + 4^ 



X 1- fix 
1 + (1-X) 2 



X 



X 



1 — x p 2 
^ x 2x 



- 4 & + ^ 



G ( «\x,p)=C F 



X X 

„2 



P- 



p + x 2 + 2x - 4 1 + /3. 



In 



+ 2/3,(2 + p) 



1 — rr 



(40) 



In order to establish the subtraction terms to be included in the GM result to ensure 
matching with the ZM result in the massless limit, we need to take the limit m — > 
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in the virtual-soft and real correction terms listed above. In this limit, the vector and 
axial-vector parts coincide. Specifically, we have 



lim 

p-^0 L 



(«), 



C F ( In 2 - + In - - 2 + vr 2 

P P 



]im B!£ ) (x,p) + Rt ) (x,p) 
p^o (1 - x)+ 



C F <J 5(1 - ar) 
1 



. 2 4 14 1 7T 

In 2 - + - In 



1 — x 
- (l + x 2 



p 2 p 2 3 
4 x 2 
P 2 
1 + x 2 



(l + x 2 ) In - -4x + ^- 
ln(l — x 



1 — x 



1 — x 



lnx 



lim 

p^O 



2C, 



1 + (1 -x) 



X 



G$\x,p)+G$\x,p) 
Comparison of Eqs. (iTTj) and (jHJ) with Eq. (jSJ) yields 

4 u) (x,p) + 4 u) (x,p) 



4 



Hm^(l-x) 4" ) (p)+5T(p) 

p^O 



+ 



(1-x) 



In - + ln(l - x) - 1 

. P 



P/ 



(41) 



x) 



lim 

p-^0 



G?(x,p) + GT(x,p) 



(«), 



where H"8lfT3 



< W/) = P^x) In ^ + C F (1 + x 2 ) <j 5(1 - x) 



dg\x,fif) = 2P$!$(x) ln^ -2C F (21nx+ 1). 



+4 i )(x,p / ), 

p f 

+ 4 1 )(x,p / ), (42) 

_J_\ 2 r in(i-x) - 

1 — xy + 1 — x + 

(43) 
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Figure 1: Kinematically accessible region in the (x, y) plain for e + e~ — > b + X and 
subsequent decay 6 — > if c + X with collinear emission of H c assuming = mz, m = 
5.0 GeV, and m H = 2.01 GeV (solid line). For comparison, also the region used in our 
analysis is shown (dashed line). 
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Figure 2: x p distributions of e + e~ — > D + + X at = 10.52 GeV evaluated in the GM 
approach with the FFs from the joint fit to the Belle [9] and CLEO pH] data including 
ISR corrections (solid line) and corresponding result with the latter subtracted (dashed 
line). 
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Figure 3: Normalized x distributions of (a) D°, (b) D + , and (c) D* + mesons from OPAL 
PI3"] and ALEPH [5 J compared to the respective fits in the GM approach from Tables HH31 
The dotted, dashed, and solid lines refer to the c-quark-initiated, 6-quark-initiated, and 
total contributions, respectively. 
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Figure 4: x p distributions of (a) D°, (b) D + , and (c) D* + mesons from Belle [9] and 
CLEO [10] compared to the respective fits in the GM approach from Tables [TH3l 
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Figure 5: x p distributions of D° mesons from (a) Belle [9] and CLEO [TO], and (b) 
normalized x distribution of D° mesons from OPAL [2J compared to the global fit from 
Table[TJ In frame (b), the dotted, dashed, and solid lines refer to the c-quark-initiated, 
6-quark-initiated, and total contributions, respectively. 




Figure 6: x p distributions of D + mesons from (a) Belle [0] and CLEO [TU], and (b) 
normalized x distribution of D + mesons from OPAL [2] compared to the global fit from 
Table [2j In frame (b), the dotted, dashed, and solid lines refer to the c-quark-initiated, 
6-quark-initiated, and total contributions, respectively. 
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Figure 7: x p distributions of D* + mesons from (a) Belle [H] and CLEO [TU], and (b) 
normalized x distributions of D* + mesons from ALEPH [3], and OPAL [3] compared to 
the global fit from Table El In frame (b), the dotted, dashed, and solid lines refer to the 
c-quark-initiated, 6-quark-initiated, and total contributions, respectively. 




Figure 8: (a) c -> D+ and (b) g -> D + FFs of the Belle/CLEO fits at [L S = 10.52 GeV as 
functions of z in the GM approach (solid lines) and in the ZM approaches with m# ^ 
(dashed lines) and m# = (dotted lines). 
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Figure 9: x p distributions of e + e — > -D + + X at y's = 10.52 GeV in the GM approach 
(solid line) and in the ZM approaches with mn 7^ (dashed lines) and m# = (dotted 
lines), all evaluated with the FFs from the Belle/CLEO-GM fit. 




Figure 10: (a) c — ► D + and (b) g — > D + FFs of the global fit in the GM approach at 
fMf — 10.52 GeV (solid lines) compared with their counterparts from Refs. pQ (dashed 
lines) and [6] (dotted lines). 
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